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This  report  describes  the  results  of  a  theoretical  investigation  into  the 
nonlinear  optical  response  of  granular  composites.  In  particular,  we  have 
investigated  the  degree  to  which  the  cubic  nonlinear  response  can  be  enhanced 
in  such  composites,  relative  to  the  same  response  in  pure  materials.  Cubic 
nonlinear  response  is  potentially  of  great  practical  importance  in  developing 
novel  materials  with  a  strongly  intensity-dependent  transmission  coefficient, 
index  of  refraction,  etc.  We  find  that  there  are  many  types  of  composite 
geometries  which  give  rise  to  such  strong  enhancements.  These  include  composites 
near  the  percolation  threshold,  composites  in  which  one  of  the  components  is 
arranged  in  the  form  of  a  fractal  cluster,  and  nonlinear  coatings  (such  as 
dyes)  adsorbed  onto  the  surface  of  particles  which  have  a  resonant  optical 
response  at  an  appropriate  frequency.  We  have  carried  out  this  investigation 
using  both  analytical  approximations  and  a  variety  of  approximate  numerical 
techniques,  including  large-scale  computer  simulations.  Finally,  we  briefly 
describe  some  work  in  progress,  in  which  a  composite  composed  of  materials  with 
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STATEMENT  OF  THE  PROBLEM 


ARO  Grant  DAAL03-92-G-0263  ran  from  February  1.  1993  to  January  31.  1994.  The 
central  goal  of  this  project  was  to  carry  out  a  theoretical  investigation  of  ways  to  increase  the 
cubic  nonlinear  susceptibility  of  composite  materials.  A  large  cubic  nonlinearitv  is  equivalent 
to  a  material  with  a  strongly  intensity-dependent  dielectric  function.  Such  materials  have 
obvious  practical  potential,  as  intensity-dependent  filters,  for  example.  In  addition,  they  are 
of  basic  interest  as  materials  with  unusual  nonlinear  optical  properties. 

It  has  recently  been  discovered  that  cubic  nonlinearities  of  this  kind  are  greatly  enhanced 
in  composite  materials.  The  central  reason  for  this  enhancement  is  that  the  local  electric  field 
can  be  greatly  enhanced  in  a  composite  material.  Hence,  a  susceptibility  which  depends  on 
a  higher  power  of  the  electric  field,  such  as  a  cubic  nonlinear  susceptibility,  can  be  increased 
even  more  than  a  linear  susceptibility.  The  goal  of  this  project  was  to  investigate  various  ways 
of  increasing  this  enhancement.  Of  particular  interest  was  the  nonlinear  response  of  certain 
dyes  coated  on  a  substrate  with  a  resonant  surface  plasmon  excitation.  In  particular,  we 
wanted  to  show  that  this  response  is  greatly  enhanced  at  frequencies  close  to  this  excitation 
(this  expectation  has  been  realized;  see  below). 


SUMMARY  OF  MOST  IMPORTANT  RESULTS 

Our  most  important  results  are  summarized  in  the  publications  which  emerged  from 
this  grant.  For  convenience,  we  also  summarize  these  results  briefly  in  the  next  several 
paragraphs,  which  consist  of  abstracts  of  the  publications  generated  with  support  from  the 
ARO  grant. 

1.  “Cubic  Nonlinearities  in  Small-Particle  Composites:  Local-Field  Induced  Giant  En¬ 
hancements,”  by  D.  Stroud  and  X.  Zhang  (Physica  A207,  pp.  55-64  (1994)).  In 
this  paper,  we  discuss  methods  for  enhancing  the  cubic  nonlinear  susceptibility  Xe 
of  a  composite  material,  starting  from  an  exact  relation  between  Xe  and  the  fourth 
moment  of  the  electric  field  in  the  related  linear  composite.  At  zero  frequencies,  in 
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a  random  metal-insulator  composite.  \e  can  be  greatly  enhanced  near  a  percolation 
threshold.  Similar  enhancement  is  obtained  when  a  linear  fractal  is  embedded  in  a 
nonlinear  host.  In  a  random  Drude  metal/insulator  composite  x,  is  greatly  enhanced 
near  surface-plasmon  resonances;  the  enhancement  shows  very  strong  structure  which 
is  nearly  undetectable  in  the  linear  response.  In  a  suspension  of  spheres  coated  with  a 
nonlinear  material,  Xe/Xcoat  >  1  near  the  surface-plasmon  resonance  frequency  of  the 
core  particle. 

2.  “Numerical  Studies  of  the  Nonlinear  Properties  of  Composites,”  by  X.  Zhang  and 
D.  Stroud  (Phys.  Rev.  B49,  944  (1994)).  In  this  paper,  we  use  both  numerical 
and  analytical  techniques  to  investigate  various  ways  of  enhancing  the  cubic  non¬ 
linear  susceptibility  Xe  of  a  composite  material.  We  start  from  the  exact  relation 
Xe  —  ZiPiXi  <  {E  ■  E)2  >i,Un  / Eq ,  where  Xi  and  Pi  are  the  cubic  nonlinear  suscep¬ 
tibility  and  volume  fraction  of  the  ith  component,  E0  is  the  applied  electric  field,  and 
<  Ea  >itiin  is  the  expectation  value  of  the  electric  field  in  the  ith  component,  calculated 
in  the  linear  limit,  where  Xi  —  0-  In  our  numerical  work,  we  represent  the  composite  by 
a  random  resistor  or  impedance  network,  calculating  the  electric  field  distributions  by 
a  generalized  transfer-matrix  algorithm.  Under  certain  conditions,  we  find  that  Xc  is 
greatly  enhanced  near  the  percolation  threshold.  We  also  find  a  large  enhancement  for 
a  linear  fractal  in  a  nonlinear  host.  In  a  random  Drude  metal/insulator  composite  %e  is 
hugely  enhanced  especially  near  frequencies  which  correspond  to  the  surface-plasmon 
resonance  spectrum  of  the  composite.  At  zero  frequency,  the  random  composite  re¬ 
sults  are  reasonably  well  described  by  a  nonlinear  effective- medium  approximation. 
The  finite-frequency  enhancement  shows  very  strong  reproducible  structure  which  is 
nearly  undetectable  in  the  linear  response  of  the  composite,  and  which  may  possibly 
be  described  by  a  generalized  nonlinear  effective- medium  approximation.  The  fractal 
results  agree  qualitatively  with  a  new  nonlinear  differential  effective- medium  approxi¬ 
mation.  Finally,  we  consider  a  suspension  of  coated  spheres  embedded  in  a  host.  If  the 
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coating  is  nonlinear,  we  show  that  \£  Xcoat  1  near  the  surface-plasmon  resonance 
frequency  of  the  core  particle. 

3.  “Effective  Linear  and  Nonlinear  Response  of  Fracta.  Clusters,  by  P.  M.  Hui  and  D. 
Stroud  (Phys.  Rev.  B49.  11729  (1994))/  In  this  paper,  we  generalize  a  differential  ef¬ 
fective  medium  approximation  previously  developed  tor  tne  effective  linear  response  of 
a  fractal  cluster,  to  treat  the  effective  nonlinear  response.  We  consider  a  material  hav¬ 
ing  a  nonlinear  J-E  characteristic  which  forms  a  cluster  in  another,  possibly  nonlinear 
material.  A  differential  equation  for  the  cluster-size-dependent  nonlinear  susceptibility 
is  derived  and  a  general  solution  is  obtained.  Numerical  examples  are  given  for  two 
special  cases:  a  cluster  of  a  nonlinear  material  in  a  linear  host,  and  a  cluster  of  a 
linear  material  in  a  nonlinear  host.  An  enhancement  in  the  nonlinear  response  due  to 
clustering  is  found  for  the  case  of  clustering  of  a  nonlinear  material  in  a  linear  host. 

4.  “Theory  of  Optical  Bistability  in  a  Weakly  Nonlinear  Composite  Medium,”  by  David 
J.  Bergman,  Ohad  Levy,  and  D.  Stroud  (Phys.  Rev.  B49,  129  (1994).  This  paper 
presents  a  new  variational  approach  to  discuss  bistable  behavior  in  composite  media 
where  one  component  is  a  nonlinear  dielectric  while  the  other  has  a  field  indepen¬ 
dent  complex  dielectric  constant  with  a  negative  real  part.  Under  certain  conditions, 
the  physical  parameters  may  be  adjusted  so  that  bistability  appears  even  though  the 
nonlinear  behavior  is  everywhere  weak,  thus  allowing  it  to  be  treated  as  a  small  per¬ 
turbation  to  the  leading  linear  behavior. 

5.  “Effective  Dielectric  Response  of  Nonlinear  Composites.'  by  K.  W.  Yu,  P.  M.  Hui,  and 
D.  Stroud  (Phys.  Rev.  B47,  14150  (1993)).  In  this  paper,  a  perturbative  approach  is 
developed  to  compute  the  local  field  for  the  case  of  a  nonlinear  inclusion  embedded  in 
a  linear  host.  The  result  is  applied  to  nonlinear  composites.  General  formulas  for  cal¬ 
culating  the  effective  nonlinear  susceptibility  up  to  the  case  of  fifth-order  nonlinearity 
are  given.  The  formulation  is  applied  to  problems  in  two  and  three  dimensions  (2D 
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and  3D).  For  2D  problems,  we  study  cylindrical  and  concentric  cylindrical  inclusions. 
By  means  of  an  exact  mapping,  the  concentric  cylinder  problem  can  mapped  onto 
that  of  an  elliptic  cylinder.  A  general  expression  is  derived  for  the  effective  nonlinear 
susceptibility  of  a  dilute  composite  of  randomly  oriented  elliptic  cylinders  embedded 
in  a  linear  host.  For  3D  problems,  we  study  the  cases  of  spherical  and  coated  spherical 
inclusions.  General  expressions  for  the  effective  nonlinear  susceptibility  are  derived  in 
the  dilute  limit  up  to  the  case  of  fifth  order  nonlinearity.  For  composites  consisting  of 
spherical  inclusions  coated  by  a  nonlinear  material  and  embedded  in  a  linear  host,  it 
is  possible  to  enhance  the  nonlinear  response  of  the  composite  by  tuning  material  pa¬ 
rameters  such  as  the  linear  dielectric  constants  of  the  host,  coating,  and  core  materials, 
and  by  adjusting  the  thickness  of  the  coating. 

6.  “Scaling  Behavior  and  Surface- Plasmon  Resonances  in  a  Model  Three-Dimensional 
Metal- Insulator  Composite,”  by  X.  Zhang  and  D.  Stroud  (Phys.  Rev.  B48,  6658 
(1993)).  In  this  paper,  we  calculate  the  a.  c.  dielectric  function  of  a  model  Drude 
metal/insulator  composite,  using  a  three-dimensional  transfer-matrix  algorithm.  The 
real  part  of  the  effective  conductivity,  Re  ge(w),  reveals  (i)  a  Drude  peak  which  appears 
only  above  the  percolation  threshold  pc;  and  (ii)  a  broad  spectrum  of  surface  plasmon 
resonances  whose  lower  edge  approaches  zero  frequency  at  pc.  Sufficiently  near  pc,  the 
dielectric  function  is  consistent  with  an  expected  scaling  form  previously  verified  in 
two  dimensions.  The  surface  plasmon  spectrum  resembles  effective-medium  predictions 
except  for  a  weak  but  persistent  peak  near  0Auj/ujp. 


WORK  IN  PROGRESS. 

Even  though  funding  has  expired,  we  are  continuing  work  along  the  lines  described  above. 
A  number  of  new  and  fascinating  lines  of  research  have  opened  up  for  possible  future  study. 
We  summarize  two  of  these  below. 
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Second-Harmonic  Generation  and  Induced  Cubic  Xonlinearities  in  Composite  Media. 
This  is  a  joint  project  with  Prof.  David  J.  Bergman  and  Ohad  Levy  of  Tel  Aviv  Uni¬ 
versity.  We  have  shown  that,  even  if  none  of  the  components  of  the  composite  have  a  cubic 
nonlinearity,  such  a  cubic  nonlinearity  can  be  indirectly  produced  m  a  medium  with  second- 
order  nonlinearity  (i.  e.  second-harmonic-generation,  or  a  nonzero  x^)-  This  induced  cubic 
nonlinearity  can  be  even  larger  than  the  enhanced  cubic  nonlinearities  analyzed  in  most  of 
the  present  work.  Thus,  this  mechanism  would  be  an  additional  means  of  producing  am 
intensity-dependent  index  of  refraction  in  a  composite  material. 

Induced  Cubic  Nonlinearities  in  a  Linear  Material.  An  even  more  intriguing  possibility 
is  that  cubic  nonlinearities  can  be  generated  in  a  material  which,  in  bulk  form,  is  strictly 
linear,  just  due  to  the  finite  size  of  the  particles  involved.  We  have  just  begun  to  investigate 
this  possibility,  which  arises  from  the  failure  of  the  quasistatic  approximation  in  composites. 
If  opportunity  arises,  we  will  continue  this  work  in  the  coming  year. 
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